AD 


Award  Number:  W81XWH-1 1-1-0127 


TITLE:  Chemical  strategy  to  translate  genetic/epigenetic  mechanisms  to  breast  cancer 
therapeutics 


PRINCIPAL  INVESTIGATOR:  Sheng  Ding,  PhD 


CONTRACTING  ORGANIZATION:  The  J.  David  Gladstone  Institutes 

Irvine,  CA  92506 


REPORT  DATE:  July  2012 


TYPE  OF  REPORT:  Annual 

PREPARED  FOR:  U.S.  Army  Medical  Research  and  Materiel  Command 
Fort  Detrick,  MD  21702-5012 


DISTRIBUTION  STATEMENT:  Approved  for  Public  Release;  Distribution  Unlimited 


The  views,  opinions  and/or  findings  contained  in  this  report  are  those  of  the  author(s)  and 
should  not  be  construed  as  an  official  Department  of  the  Army  position,  policy  or  decision  unless 
so  designated  by  other  documentation. 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and  maintaining  the 
data  needed,  and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggestions  for  reducing 
this  burden  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202- 
4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently 
valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 

1 .  REPORT  DATE  2.  REPORT  TYPE 

01-07-2012  Annual 

3.  DATES  COVERED 

July  1,2011  -  June  30,  2012 

4.  TITLE  AND  SUBTITLE 

Chemical  Strategy  to  Translate  Genetic/Epigenetic  Mechanisms  to  Breast  Cancer 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

W81XWH-1 1-1-0127 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

Sheng  Ding,  PhD 

E-Mail:  Sheng. ding® gladstone.ucsf.edu 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

The  J.  David  Gladstone  Institutes 

Irvine,  CA  92506 

8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 

9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

U.S.  Army  Medical  Research  and  Materiel  Command 

Fort  Detrick,  Maryland  21702-5012 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Approved  for  Public  Release;  Distribution  Unlimited 


13.  SUPPLEMENTARY  NOTES 


14.  ABSTRACT 

This  project  is  designed  to  use  a  molecular  signature  strategy  to  develop  new  therapeutics  against  metastatic  breast 
cancer.  We  have  recently  published  this  novel  pathway-centric  technology  developed  on  a  prostate  cancer  model  (Li  et  al., 
2012b).  The  goal  for  this  project  is  to  define  a  panel  of  genes  that  are  tightly  associated  with  breast  cancer  transition  from 
initial  epithelial  morphology  to  mesenchymal-like  cells.  This  process,  known  as  the  EMT  transition,  has  been  shown  to  play  a 
critical  role  for  breast  cancer  metastasis.  Once  the  panel  of  EMT-linked  is  defined  on  a  triple  negative  breast  cancer,  the  next 
step  is  to  conduct  small  molecular  screens  against  such  collection  of  responsive  genes.  Resulting  candidate  hits  will  be 
characterized  on  cell  and  animal  models  to  determine  their  biological  effects  in  inhibiting  breast  cancer  metastasis  to  distal 
organs.  This  approach  is  uniquely  suited  to  attack  disease  pathways  in  the  absence  of  druggable  targets  that  can  be  pursued 
by  conventional  approaches. 


15.  SUBJECT  TERMS 

No  subject  terms  provided. 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

OF  ABSTRACT 

18.  NUMBER 

OF  PAGES 

19a.  NAME  OF  RESPONSIBLE  PERSON 

USAMRMC 

a.  REPORT 

u 

b.  ABSTRACT 

U 

c.  THIS  PAGE 

U 

uu 

12 

19b.  TELEPHONE  NUMBER  (include  area 
code) 

Table  of  Contents 


Page 


Introduction .  1 

Body .  1  -  2 

Key  Research  Accomplishments .  2 

Reportable  Outcomes .  3 

Conclusion .  3 

References .  3 

Appendices .  3-9 


Introduction 


This  project  is  designed  to  use  a  m  olecular  signature  strategy  to  develop  new 
therapeutics  against  metastatic  breast  cancer.  We  have  recently  published  this  novel  pathway¬ 
centric  technology  developed  on  a  prostate  cancer  model  (Li  et  al.,  2012b).  The  goal  fro  this 
project  is  to  define  a  panel  of  genes  that  are  tightly  associated  with  breast  cancer  transition  from 
initial  epithelial  morphology  to  mesenchymal-like  cells.  This  process,  known  as  the  EMT 
transition,  has  been  shown  to  play  a  critical  role  for  breast  cancer  metastasis.  Once  the  panel  of 
EMT-linked  is  defined  on  a  triple  negative  breast  cancer,  the  next  step  is  to  conduct  small 
molecular  screens  against  such  collection  of  responsive  genes.  Resulting  candidate  hits  will  be 
characterized  on  cell  and  animal  models  to  determine  their  biological  effects  in  inhibiting  breast 
cancer  metastasis  to  distal  organs.  This  approach  is  uniquely  suited  to  attack  disease  pathways 
in  the  absence  of  druggable  targets  that  can  be  pursued  by  conventional  approaches. 


Body  of  the  Progress  Report 

According  to  our  original  Statement  of  Work,  the  project  is  to  proceed  in  three  phases. 

Aim  1:  Define  the  gene  signature  associated  with  the  EMT  program 

(1)  Perform  RNAi  on  SUM1315  cells  against  Twist,  Snail,  Slug,  and  SIP1 

(2)  Confirm  the  RNAi  efficiency  by  Western  blotting 

(3)  Conduct  RNA-seq  analysis  on  both  wt  and  RNAi  knockdown  cells 

(4)  Conduct  bioinformatics  analysis  of  the  RNA-seq  data  to  define  regulated  genes  in  each  case 

(5)  Deduce  specific  and  common  gene  signature  trigger  by  individual  EMT  inducers 

(6)  Validate  a  subset  of  regulated  genes  by  real  time  PCR 

Aim  1  has  been  completed  entirely  in  our  partner  lab,  Dr.  XiangDong  Fu,  at  UCSD. 


Aim  2:  Screen  for  small  molecules  that  can  specifically  repress  the  EMT  program 

In  close  collaboration  with  Dr.  Fu  lab  at  UCSF,  we  participated  in  developing  the  RASL- 
seq  platform  for  screening  small  molecules  in  breast  cancer  cells  that  can  suppress  the  EMT 
gene  signature.  From  expression  analysis,  205  genes  were  selected  (202  common  genes  + 
Twistl,  Sipl  and  Slug)  to  design  RASL-seq  oligo  pairs,  and  subsequently  experimentally 
validated  to  (1)  eliminate  those  that  give  rise  to  too  many  counts  (this  will  consume  our 
sequence  space),  (2)  compare  the  results  with  our  genome-wide  and  RT-qPCR  results  to 
ensure  that  individual  oligo  pairs  report  the  anticipated  effects  in  regulated  gene  expression,  (3) 
re-design  and  re-synthesize  those  that  are  not  efficient  in  targeting  intended  genes,  and  (4) 
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conduct  initial  screen  on  our  robot  to  determine  if  the  cell  culture  system  and  the  oligo  pool  can 
be  adapted  to  high  throughput  operation. 


In  addition,  we  spent  a  major  effort  on  optimizing  the  cell  culture  and  treatment 
conditions  for  conducting  the  RASL-seq  screen.  The  optimization  in  384-well  format  included 
examining  plating  cell  density,  media  (with  serum  and  additional  supplemented  growth  factors), 
drug  treatment  concentration,  treatment  time,  media  change  etc,  so  that  cells  would  have  the 
best  response  and  samples  could  be  processed  in  parallel  for  RASL-seq.  After  those 
optimizations,  we  carried  out  the  screening  with  4000  compounds.  The  initial  data  is  currently 
being  analyzed.  If  things  go  well,  we  will  start  full  deck  screening  effort  in  the  following  months 
and  start  doing  more  chemistry  work  and  mechanistic  studies. 


Aim  3:  Test  the  effect  of  drug  leads  on  cell  and  animal  models 

(1)  Conduct  Western  blotting  of  epithelial  and  mesenchymal  markers  on  SUM1315,  Bt594,  and 
MDA157  cells  before  and  after  treatment  with  drug  leads 

(2)  Perform  cell  migration  and  invasion  assay  on  mock-treated  and  drug-treated  cells 

(3)  Test  the  effect  of  leading  candidates  on  breast  cancer  metastasis  on  nude  mice 


This  aim  will  be  pursued  in  the  second  year  once  screening  and  follow-up  confirmation  are 
completed 


Key  Research  Accomplishments 

•  Conducted  RNAi  on  a  breast  cancer  model 

•  Performed  genome-wide  analysis  to  define  the  EMT-linked  molecular  signature 

•  Designed,  synthesized,  and  tested  initial  oligo  pool  for  chemical  screening 

•  Tested  the  screening  procedure  on  our  assay  robot 

•  Adapted  the  cell  model  on  screening  robot 

•  Tested  and  optimized  conditions  for  chemical  screening,  and  carried  out  initial  test  run 
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Reportable  Outcomes 

The  screening  technology  platform  paper  was  published  in  PNAS  (listed  below).  In  addition,  we 
have  been  presenting  our  technology  development  efforts  in  a  recent  Society  for  Laboratory 
Automation  and  Screening  (SLAS)  meeting  in  San  Diego. 


Conclusion 

We  believe  that  we  are  progressing  exactly  as  we  described  in  our  proposal.  We  encountered 
some  problems  in  validating  some  of  initially  synthesized  oligos,  which  might  be  related  to  pool 
quality  of  chemical  synthesis.  We  are  now  addressing  these  problems.  Other  than  this  technical 
issue,  we  are  on  target  to  our  set  goals. 


References 

None 
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The  advent  of  powerful  genomics  technologies  has  uncovered 
many  fundamental  aspects  of  biology,  including  the  mechanisms  of 
cancer;  however,  it  has  not  been  appropriately  matched  by  the 
development  of  global  approaches  to  discover  new  medicines 
against  human  diseases.  Here  we  describe  a  unique  high-through- 
put  screening  strategy  by  high-throughput  sequencing,  referred  to 
as  HTS2,  to  meet  this  challenge.  This  technology  enables  large-scale 
and  quantitative  analysis  of  gene  matrices  associated  with  specific 
disease  phenotypes,  therefore  allowing  screening  for  small  mole¬ 
cules  that  can  specifically  intervene  with  disease-linked  gene-ex¬ 
pression  events.  By  initially  applying  this  multitarget  strategy  to 
the  pressing  problem  of  hormone-refractory  prostate  cancer,  which 
tends  to  be  accelerated  by  the  current  antiandrogen  therapy,  we 
identify  Peruvoside,  a  cardiac  glycoside,  which  can  potently  inhibit 
both  androgen-sensitive  and  -resistant  prostate  cancer  cells  with¬ 
out  triggering  severe  cytotoxicity.  We  further  show  that,  despite 
transcriptional  reprogramming  in  prostate  cancer  cells  at  different 
disease  stages,  the  compound  can  effectively  block  androgen  re¬ 
ceptor-dependent  gene  expression  by  inducing  rapid  androgen  re¬ 
ceptor  degradation  via  the  proteasome  pathway.  These  findings 
establish  a  genomics-based  phenotypic  screening  approach  capa¬ 
ble  of  quickly  connecting  pathways  of  phenotypic  response  to  the 
molecular  mechanism  of  drug  action,  thus  offering  a  unique  path¬ 
way-centric  strategy  for  drug  discovery. 

chemical  screening  |  gene  signature 

It  is  of  utmost  importance  to  match  the  power  of  functional 
genomics  in  interrogating  diseased  cells/tissues  with  potent 
drug-discovery  approaches.  Although  target-centric  approaches 
have  been  favored  in  the  past  decade,  phenotypic  screening 
appears  to  have  out-paced  such  mechanism-based  screening 
strategies  in  discovering  “first-in-class”  drugs,  thus  igniting  recent 
debate  on  the  merit  of  target-based  strategies  (1).  This  debate  is 
important  because  analysis  of  US  Food  and  Drug  Administra¬ 
tion-approved  drugs  in  recent  decades  have  revealed  low  pro¬ 
ductivity  in  drug  research  and  development,  despite  staggering 
investment  in  the  pharmaceutical  industry  (2). 

Although  phenotypic  approaches  score  the  final  functional 
outcomes,  it  is  challenging  to  optimize  candidate  drugs  without 
knowing  their  mechanism  of  action  and  many  procedures  have 
limited  capacity  in  implementing  high-throughput  screening.  In 
contrast,  target-centric  approaches  have  their  own  problems  be¬ 
cause  specific  molecular  hypotheses  based  on  the  existing 
knowledge  may  or  may  not  be  related  to  disease  phenotype.  A 
proposed  solution  to  these  problems  is  to  monitor  the  collective 
response  of  all  relevant  genes  to  a  specific  disease  phenotype  (3), 
but  this  has  been  a  major  challenge  with  any  existing  technologies. 

A  “quick-win/fast-fail”  strategy  has  been  proposed  to  streamline 
initial  candidate  hits  in  early  phases  to  offset  high  attrition  rates  in 
drug  discovery  (4).  This  strategy  begs  the  question  of  how  to  retain 
the  advantages  of  both  phenotypic  and  target-centric  screening 
approaches  to  quickly  lead  to  promising  drug  candidates.  One  idea 
is  to  take  molecular  approaches  to  conduct  “phenotypic  screen¬ 
ings”  by  using  a  set  of  genes  to  report  a  disease  phenotype,  thus 
allowing  screening  for  small  molecules  that  can  effectively  block 


disease-associated  gene-expression  events.  Several  proof-of-con- 
cept  experiments  for  this  gene-signature  approach  have  been 
performed  on  a  number  of  disease  paradigms  by  PCR  (5)  or  hy¬ 
bridization  on  Luminex  beads  (6,  7).  However,  none  of  these 
methods  has  yet  reached  the  throughput  or  cost-effectiveness  that 
is  required  for  large-scale  applications  in  drug  discovery. 

Here,  we  report  the  development  of  a  pathway-centric  high- 
throughput  screening  strategy  by  taking  advantage  of  the  ever- 
increasing  power  of  high-throughput  sequencing  (this  strategy  is 
hereafter  referred  to  as  HTS2).  The  technology  permits  quanti¬ 
tative  analysis  of  a  gene  matrix  directly  in  cell  lysates.  We  have 
initially  applied  this  technology  on  a  prostate  cancer  cell 
(LNCaP)  model  to  identifying  small  molecules  that  can  block  the 
androgen  receptor  (AR)-mediated  gene  expression,  because  AR 
overexpression  has  been  linked  to  prostate  cancer  progression  to 
androgen-refractory,  incurable  tumors  (8,  9).  By  following  the 
AR  pathway,  we  identified  several  classes  of  compounds,  one  of 
which  belongs  to  cardiac  glycosides  currently  used  for  treating 
congestive  heart  failure  and  arrhythmias,  but  also  known  for 
their  broad  anticancer  activities  on  cancer-cell  models  (10,  11). 
Interestingly,  a  recent  epidemiological  study  revealed  that  long¬ 
term  use  of  cardiac  glycosides  has  significant  benefit  in  pre¬ 
venting  prostate  cancer  (12).  We  now  showed  that  this  activity  is 
likely  a  result  of  the  ability  of  cardiac  glycosides  to  cause  AR 
destabilization,  thereby  effectively  blocking  AR-dependent  gene 
expression  and  cell  proliferation  on  both  androgen-sensitive  and 
-resistant  prostate  cancer  cells.  These  findings  validate  the  po¬ 
tential  of  this  recently  developed  HTS2  technology  in  pathway¬ 
centric  chemical  screenings,  matching  the  advances  in  functional 
genomics  to  the  development  of  new  anticancer  therapies. 

Results 

HTS  Platform  Based  on  Next-Generation  Sequencing.  As  diagrammed 
in  Fig.  L4,  we  first  profile  gene  expression  in  a  chosen  cell  type  to 
define  a  panel  of  genes  associated  with  a  disease  phenotype.  To 
quantify  these  genes  in  a  high-throughput  manner,  we  use  the 
RNA  annealing,  selection,  ligation  (RASL)  strategy,  originally 
designed  to  profile  mRNA  isoforms  using  pooled  pairs  of  oligo¬ 
nucleotides,  each  flanked  by  a  universal  primer  to  target  specific 
splice  junctions  in  spliced  mRNAs  (13,  14).  Upon  annealing  to 
total  RNA  followed  by  solid-phase  selection  via  oligo-dT  or  bio¬ 
tinylated  total  RNA,  paired  DNA  probes  templated  by  specific 
RNA  sequences  can  be  ligated  by  T4  DNA  ligase,  thus  converting 
singleton  probes  to  PCR  amplicons.  Minimal  bias  is  introduced  in 
this  process  because  of  the  uniform  length  and  relatively  balanced 


Author  contributions:  H.L.,  M.G.R.,  S.D.,  and  X.-D.F.  designed  research;  H.L.,  H.Z.,  D.W., 
and  J.Q.  performed  research;  H.L.  and  X.L.  contributed  new  reagents/analytic  tools;  H.L., 
D.W.,  Y.Z.,  and  X.-D.F.  analyzed  data;  and  H.L.,  M.G.R.,  S.D.,  and  X.-D.F.  wrote  the  paper. 

The  authors  declare  no  conflict  of  interest. 

Data  deposition:  The  RNA-seq  data  reported  in  this  paper  have  been  deposited  in  the  Gene 
Expression  Omnibus  (GEO)  database,  www.ncbi.nlm.nih.gov/geo  (accession  no.  GSE35126). 

To  whom  correspondence  may  be  addressed.  E-mail:  mrosenfeld@ucsd.edu,  sheng.ding@ 
gladstone.ucsf.edu,  or  xdfu@ucsd.edu. 

This  article  contains  supporting  information  online  at  www.pnas.org/lookup/suppl/doi:  10. 
1073/pnas.1200305109/-/DCSupplemental. 


www.pnas.org/cgi/doi/1 0. 1 073/pnas.  1 200305109 


4 


PNAS  |  March  20, 2012  |  vol.  109  |  no.  12  |  4609-4614 


MEDICAL  SCIENCES 


'  Gene  1 
Gene  2 


Solid  phase  selection 
followed  by  oligo 
ligation  on  RNA  template 


Amplify  ligated  products 
with  bar-coded  primers 

Bar-code 


(A)n 

(A)n 


(A)n 


B 


DHU 


DMSQ  DHT  Cyp  Flu  CDX 


Average  ES: 


0.24  0.65 


±0.01  ±0.06 
log2  -4  -2  0  +2  +4 


0.64 

±0.06 


C 

DMSQ  DHT  Cyp 


Multiplex  Sequencing,  first  in  the  ligated 
oligo  region  and  then  in  the  bar-code  regions 


Average  SS: 


0.98  0.93  0.92 
±0.04  ±0.04  ±0.07 


Fig.  1.  The  scheme  of  the  HTS2  technology.  (>4)  The  flow  of  the  HTS2 
technology.  ( B  and  C)  Representative  responses  to  androgen  (DHT)  and  the 
effects  of  known  antiandrogen  compounds  on  a  panel  of  DHT-responsive 
genes  ( B )  and  a  set  of  housekeeping  and  cytotoxicity-related  genes  (C). 
Three  vehicle-treated  samples  are  averaged  to  serve  as  the  baseline  and 
variations  from  the  baseline  are  color-coded:  red  for  DHT-induced  genes  and 
blue  for  DHT-suppressed  genes.  The  ES,  SS,  and  SD  are  shown  on  the  bottom 
of  each  panel  for  each  compound. 


GC  content  of  designed  probes  (15).  Instead  of  quantifying  the 
amplified  products  on  microarrays,  as  before,  we  can  now  use 
deep-sequencing  to  directly  count  the  correctly  ligated  products. 
By  using  bar-coded  primers  (each  contains  a  unique  7-nt  se¬ 
quence),  multiple  samples  can  be  pooled  for  parallel  quantifica¬ 
tion  by  first  reading  target  sequences  followed  by  sequencing  the 
bar-coded  region  to  decode  different  samples  in  the  pool. 

This  experimental  scheme  is  fully  amenable  to  automation  and 
direct  transcript  analysis  in  the  cell  lysate,  two  critical  parameters 
for  high-throughput  applications.  As  hybridization  can  take  place 
in  the  presence  of  detergent  and  high  salt,  the  annealing  step  is 
fully  compatible  to  standard  hybridization  conditions.  Biotinylated 
oligo-dT  included  in  the  annealing  step  captures  spliced  mRNA 
from  the  cell  lysate  along  with  annealed  probes  on  them  via 
streptavidin  selection.  RNA  capture  can  be  alternatively  per¬ 
formed  on  oligo-dT-coated  plates,  which  produced  similar  results. 
All  subsequent  washing  and  ligation  steps  are  carried  out  on  the 
solid  phase.  Ligated  products  are  released  in  H20,  converted  to 
bar-coded  amplicons  by  PCR,  and  pooled  for  deep  sequencing.  We 
have  fully  implemented  this  HTS2  strategy  on  a  Biomek  FX  robot. 

Robust  Performance  of  HTS2  on  Androgen-Regulated  Genes.  To 

demonstrate  this  unique  technology,  we  initially  selected  two 
dozen  androgen-responsive  genes  for  comparison  between  fold- 
differences  detected  by  HTS2  and  RT-qPCR  on  mock-treated  and 
androgen  (dihydrotestosterone,  DHT)-treated  LNCaP  cells,  which 
validated  the  ability  of  the  HTS2  technology  in  detecting  quanti¬ 
tative  differences  (Fig.  SL4).  We  noticed  that  different  probe  sets 


against  the  same  genes  often  showed  differences  in  annealing/li¬ 
gation,  thus  resulting  in  different  tag  counts,  but  fold-differences 
were  relatively  consistent  (Fig.  SIB).  This  finding  actually  allowed 
us  to  select  low-efficiency  probe  sets  for  abundant  transcripts  and 
high-efficiency  probe  sets  for  low-copy  transcripts  to  balance  the 
sequence  space  for  transcripts  of  varying  abundance. 

By  requiring  at  least  200  counts  per  transcript  in  a  panel  con¬ 
sisting  of  ~100  targets,  we  initially  tested  pooling  all  samples  in 
384  wells  for  sequencing  in  a  single  lane  on  the  Illumina  GAII 
system,  which  normally  generates  ~20  million  reads,  thus  suffi¬ 
cient  for  the  required  count  capacity  (200  x  100  x  384  =  ~8 
million  reads).  More  recently,  we  further  successfully  tested  an 
even  higher  multiplexicity  by  pooling  samples  from  four  such  384- 
well  plates  for  sequencing  on  HiSeq2000,  which  has  the  capacity 
to  produce  up  to  200  million  counts  per  lane.  Thus,  on  a  pair  of 
eight-lane  flowcells  in  a  single  HiSeq2000  run,  the  system  permits 
screening  >20  thousand  samples  (16  x  4  x  384  =  24,576). 

Application  of  HTS2  to  Analysis  of  Antiandrogen  Activities.  We  next 
applied  the  HTS2  platform  to  identify  small  molecules  that  can 
effectively  block  androgen-induced  gene  expression  in  LNCaP 
cells.  Unlike  early  efforts  based  on  a  single  expression  reporter 
(16),  we  selected  a  panel  of  70  androgen-induced  (both  up-  and 
down-regulated)  genes  to  represent  the  AR  pathway,  determined 
by  RNA-seq  (17).  We  also  included  30  housekeeping  and  cyto¬ 
toxicity-related  genes  (18)  as  specificity  controls  (Table  SI).  The 
robotic  RASL  assay  demonstrated  high  reproducibility  among 
biological  repeats  of  mock- treated  and  DHT-treated  samples 
(Fig.  S2  A  and  B)  and  the  ability  to  detect  anticipated  androgen 
responses  (Fig.  S2C).  We  then  tested  several  known  androgen 
antagonists,  including  Cyproterone  (Cyp),  Flutamide  (Flu),  and 
Bicalutamide  (CDX).  Cyp  showed  a  modest  effect  but  both  Flu 
and  CDX  were  more  potent  than  Cyp  in  blocking  DHT-induced 
gene  expression  (Fig.  IB),  all  of  which  had  minimal  effect  on  the 
built-in  controls  (Fig.  1C). 

Previous  works  used  summed  gene  expression,  k-nearest 
neighbors,  or  naive  Bayes  classification  to  quantify  the  anti¬ 
androgen  effect  (5).  These  methods  generate  an  index  to  each 
compound  based  on  alteration  in  gene  expression  toward  a  desired 
direction.  Using  a  panel  of  controls  to  simultaneously  evaluate  the 
specificity,  we  generated  both  the  efficacy  and  specificity  scores  for 
each  compound  by  using  summed  gene  expression  {Methods).  The 
Efficacy  score  (ES)  measures  independent  contribution  of  each 
gene  to  antagonizing  the  DHT  effect,  which  is  summed  to  produce 
a  score  ranging  from  -1  (for  full  effect  in  enhancing  the  DHT 
effect)  to  +1  (for  full  effect  in  suppressing  the  DHT  effect).  The 
Specificity  score  (SS)  determines  the  number  of  control  genes  that 
remain  unaltered  (twofold),  ranging  from  0  (for  nonspecific  effect) 
to  +1  (for  lack  of  nonspecific  effect).  Therefore,  we  are  able  to 
generate  a  pair  of  numerical  indexes  for  each  of  the  known  anti¬ 
androgen  compounds  (Fig.  1 B  and  C).  The  data  fully  concur  with 
different  potencies  of  these  drugs  in  treating  prostate  cancer 
patients  in  the  clinic  (19),  demonstrating  the  ability  of  the  HTS2 
platform  in  characterizing  androgen  antagonists  based  on  the 
collective  response  of  androgen-responsive  genes. 

Identification  of  Unique  Antiandrogen  Compounds.  On  the  estab¬ 
lished  HTS2  platform,  we  next  screened  a  collection  of  compounds 
(~4,000),  consisting  of  human-experienced  drugs,  a  set  of  natural 
products,  and  various  synthetic  inhibitors  of  known  enzymes.  We 
included  CDX  as  a  benchmark  for  antiandrogen  activity  and  full- 
effect  mimics  (LNCaP  cells  not  treated  with  DHT  or  compound, 
thus  mimicking  complete  suppression  of  the  DHT  effect).  As 
expected,  most  full-effect  mimics  (blue  dots)  showed  high  ES  and 
SS  scores;  DHT-treated  samples  (pink  dots)  produced  low  ES  and 
high  SS  scores;  and  treatments  with  CDX  (red  dots)  generated 
relatively  high  ES  and  SS  scores  (Fig.  24  and  Fig.  S3).  In¬ 
terestingly,  although  most  test  compounds  (green  dots)  lacked  any 
antiandrogen  effects  (low  ES  scores),  we  identified  a  number  of 
compounds  that  exhibited  high  ES  and  SS  scores,  indicating  that 


5 


Li  et  al. 


4610  |  www.pnas.org/cgi/doi/1 0.1 073/pnas.l 2003051 09 


o 

CD 

O 

it 

LU 


Specificity 


Bicalutamide 


0.5nM  5nM  50nM  0.5nM  5|jM  50hM 


Flutamide 


Q  Group  1  Group  2  Group  3  Group  4 


+1 


Fig.  2.  Identification  of  unique  anti-AR  compounds.  (A)  Two-dimensional  plot  of  screened  compounds  (plots  of  individual  plates  are  shown  in  Fig.  S3).  Blue,  full- 
effect  mimics  (no  DHT  treatment  and  no  compound;  pink,  DHT  treatment  alone;  red,  CDX  on  DHT-treated  cells;  and  green,  compound  on  DHT-treated  cells.  ( B ) 
Clustering  analysis  of  top  candidate  hits.  Red  represents  the  effect  on  suppressing  DHT-induced  genes;  green  shows  the  effect  on  restoring  DHT-repressed  genes.  The 
compounds  identified  from  published  screenings  as  described  in  the  text  are  highlighted  in  green,  blue,  and  pink;  cardiac  glycosides  are  labeled  in  red.  (C)  Titration 
and  deduced  ESmax  and  IC50  for  individual  compounds.  SDs  are  based  on  triplicate  measurements.  Compound  structure  is  shown  on  top  of  each  titration  curve. 


these  candidate  hits  may  function  similarly  or  better  than  CDX  in 
inhibiting  DHT-responsive  genes  in  LNCaP  cells. 

The  anti-DHT  effects  of  top  candidates  are  displayed  by  hi¬ 
erarchical  clustering  based  on  their  impact  on  individual  genes  in 
the  panel  where  red  represents  suppression  of  DHT-induced 
gene  expression  ranging  from  0  (no  effect)  to  + 1  (full  suppres¬ 
sion),  whereas  green  indicates  restoration  of  DHT-repressed 
gene  expression,  which  ranges  from  0  (no  effect)  to  -1  (full 
restoration)  (Fig.  2 B).  This  analysis  revealed  four  groups:  groups 
1  and  4  showed  effective  suppression  of  DHT-induced  genes,  but 
little  effect  on  restoring  DHT  down-regulated  genes,  indicating 
that  these  compounds  may  impair  gene  expression  in  some 
general  ways.  In  comparison,  groups  2  and  3  antagonized  DHT- 
responsive  genes  in  both  directions. 

We  note  two  estrogens,  Estrone  and  Mestranol  (green  in  Fig. 
2 B),  in  the  top  hits,  suggesting  that  the  activation  of  the  estrogen 
pathway  might  interfere  with  the  status  of  the  androgen  pathway, 
which  is  consistent  with  the  reported  effect  of  estrogen  in 
inhibiting  prostate  cancer  cells  through  the  estrogen  receptor  p 
(ERp)  (20).  Two  compounds  Pyrvinium  Panoate  and  Exalamide 
(blue  in  Fig.  2 B)  were  previously  identified  from  an  AR  con¬ 
formational  screen  (21,  22).  We  also  identified  several  ionophors 
(pink  in  Fig.  2 B),  including  Salinomycin,  a  compound  recently 
shown  to  selectively  inhibit  breast  cancer  stem  cells  (23),  and 
Nigericin,  which  was  previously  identified  from  a  screen  based  on 
suppressing  the  prostate-specific  antigen  (PSA)  (16). 

Interestingly,  a  major  fraction  of  candidate  hits  belongs  to  the 
cardiac  glycoside  family,  including  Strophanthidin,  Ouabain, 
Proscillaridin,  Peruvoside,  Digoxin,  Concallatoxin,  and  Digitoxin 


(red  in  Fig.  2 B),  all  in  group  3.  We  selected  several  of  these 
cardiac  glycosides  to  confirm  their  anti-androgen  effects  and 
derive  the  maximum  effect  based  on  the  ES  (ESmax)  and  the  half¬ 
maximum  inhibitory  concentration  (IC50)  (Fig.  2C).  The  results 
reveal  that  cardiac  glycosides,  particularly  Peruvoside  and  Stro¬ 
phanthidin,  are  more  potent  than  CDX  and  Flu  in  blocking 
DHT-induced  gene  expression  in  LNCaP  cells  (Fig.  2C). 

Peruvodise  Potently  Inhibits  Cell  Proliferation  Without  Inducing 
Cytotoxicity.  We  next  focused  on  understanding  how  cardiac 
glycosides  suppress  AR-mediated  gene  expression  for  two  rea¬ 
sons.  First,  various  cardiac  glycosides  exhibit  broad  anticancer 
effects,  including  prostate  cancer,  both  in  vitro  (24,  25)  and  on 
castration-resistant  prostate  tumors  in  animal  models  (26).  Sec¬ 
ond,  a  recent  epidemiological  study  revealed  that  long-term  use 
of  Digoxin,  a  cardiac  glycoside  widely  used  to  treat  congestive 
heart  diseases,  significantly  reduces  the  risk  of  prostate  cancers 
(12).  Our  results  suggest  that  cardiac  glycosides  may  directly 
intervene  with  the  AR  pathway  in  prostate  cancer  cells. 

To  determine  if  the  inhibitory  effect  on  prostate  cancer  cells 
depends  on  a  functional  AR,  we  tested  several  cardiac  glycosides 
at  5  pM  (the  same  concentration  we  used  in  the  primary  screen) 
on  two  pairs  of  isogenic  prostate  cancer  cells  (Fig.  3/4).  One  pair  is 
androgen-sensitive  LNCaP  cells  and  their  androgen-resistant 
derivative  LNCaP-abl  cells;  the  other  pair  is  androgen-resistant, 
AR-negative  PC3  cells  and  AR  re-expressed  PC3-AR  cells.  All 
three  cardiac  glycosides  tested  began  to  inhibit  proliferation  of 
LNCaP  cells  within  the  first  24  h  of  treatment,  with  little  effect  on 
LNCaP-abl  cells  and  no  effect  on  PC3  and  PC3-AR  cells.  By  day  2, 
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Fig.  3.  Effects  of  lead  compounds 
on  cell  proliferation  and  apopto¬ 
sis.  (A)  Effect  of  each  compound 
(5  |iM)  on  DHT-treated  LNCaP  cells 
and  on  androgen-independent 
LNCaP-abl,  PC3,  and  PC3-AR  cells. 

(B)  Effect  of  each  compound  at 
different  concentrations  on  an¬ 
drogen-resistant  LNCaP-abl  cells. 

(C)  Apoptosis  induced  on  LNCaP 
(Left)  and  LNCaP-abl  (Right)  cells 
by  individual  compounds  with 
Nigericin  as  a  positive  control.  (D) 
Competitive  binding  determined 
with  3H-labeled  DHT  in  the  pres¬ 
ence  of  increasing  concentrations 
of  individual  cold  DHT,  CDX,  and 
Peruvoside. 


the  cardiac  glycosides  were  quite  effective  in  inhibiting  LNCaP- 
abl  cells,  but  still  lacked  effect  on  PC3  cells.  Interestingly,  PC3- 
AR  cells  appear  to  have  gained  a  degree  of  sensitivity  to  cardiac 
glycosides  (Fig.  3/4).  By  day  3,  proliferation  was  inhibited  in  all 
cell  types.  These  observations  suggest  that  a  functional  AR  may 
render  PC3  cells  sensitive  to  inhibition  by  cardiac  glycosides. 

To  further  characterize  these  cardiac  glycosides,  we  treated 
LNCaP-abl  cells  for  3  d  with  each  compound  at  different  con¬ 
centrations  (Fig.  3 B).  Peruvoside  blocked  cell  growth  at  50  nM; 
both  Digoxin  and  Strophanthidin  required  500  nM  to  achieve  the 
same  effect  on  the  androgen-resistant  LNCaP-abl  cells,  although 
Strophanthidin  effectively  inhibited  androgen-induced  gene  ex¬ 
pression  as  Peruvoside  did  on  androgen-sensitive  LNCaP  cells. 
Early  studies  suggest  that  cardiac  glycosides  inhibit  cell  pro¬ 
liferation  through  the  induction  of  apoptosis  (24-26).  We  asked 
whether  induced  apoptosis  was  sufficient  to  account  for  the  strong 
effect  of  Peruvoside  on  inhibiting  cell  proliferation.  By  monitoring 
activated  Caspases  3  and  7,  we  found  that  Peruvoside  has  a  de¬ 
tectable  degree  of  induced  apoptosis  on  LNCaP  cells,  but  no  effect 
on  LNCaP-abl  cells  (Fig.  3C).  Thus,  Peruvoside  can  effectively 
block  cell  proliferation  without  triggering  general  cytotoxic  re¬ 
sponse  on  androgen-resistant  prostate  cancer  cells.  In  contrast, 
Nigercin  induced  apoptosis  on  both  LNCaP  and  LNCaP-abl  cells. 

As  cardiac  glycosides  are  themselves  cardiotonic  steroids  (see 
their  core  steroid  structure  in  Fig.  2C),  and  in  light  of  a  recent  finding 
that  Digoxin  can  directly  bind  to  a  specific  nuclear  receptor  (RORy  t) 
in  T  cells  (27),  we  examined  whether  Peruvoside  might  compete  for 
androgen  binding.  We  found  no  evidence  for  Peruvoside  to  bind  to 
the  ligand  binding  pocket  in  AR,  because  it  could  not  compete  with 
binding  of  3H-labeled  DHT  to  endogenous  AR  on  LNCaP  cells.  In 
contrast,  unlabeled  DHT  and  CDX  could  completely  or  partially 
compete,  respectively  (Fig.  3D).  These  results  ruled  out  the  possi¬ 
bility  that  Peruvoside  acts  as  a  competitive  androgen  antagonist. 

Cardiac  Glycosides  Block  the  Entire  AR  Pathway.  Our  primary  screen 
was  conducted  on  androgen-sensitive  LNCaP  cells,  yet  these 


compounds  also  exhibited  significant  effects  on  androgen-re¬ 
sistant  LNCaP-abl  cells  that  have  a  largely  distinct  AR-dependent 
gene-expression  program  as  a  result  of  transcriptional  reprog¬ 
ramming  (28).  We  therefore  asked  how  they  might  affect  AR- 
mediated  gene  expression,  even  after  transcriptional  reprogram¬ 
ming,  by  performing  genome-wide  analysis  by  RNA-seq  (29). 

We  first  identified  AR-dependent  genes  on  LNCaP-abl  cells 
by  AR  RNAi,  and  then  compared  them  to  altered  genes  in 
cardiac  glycoside-treated  cells  (Fig.  4).  Based  on  a  stringent 
cutoff  (>twofold,  P  <  0.01),  we  identified  2,056  genes  that 
responded  to  at  least  one  treatment  (AR  knockdown  or  treat¬ 
ment  with  individual  cardiac  glycosides).  When  ranked  by  aver¬ 
aged  responses  to  cardiac  glycoside  treatments,  we  found  that 
the  induced  gene  expression  in  AR  knockdown  cells  largely 
matched  those  in  cardiac  glycoside-treated  cells  and  most  of  the 
overlapped  genes  changed  in  the  same  direction.  Taken  to¬ 
gether,  these  data  demonstrated  that  cardiac  glycosides  are  able 
to  selectively  and  effectively  block  AR-dependent  gene  expres¬ 
sion  in  LNCaP-abl  cells,  even  though  the  AR  program  has  been 
dramatically  altered  compared  with  androgen-sensitive  LNCaP 
cells  to  support  their  androgen  independent  growth. 

Peruvoside  Acts  as  the  Potent  Inducer  of  AR  Degradation.  To  further 
understand  how  cardiac  glycosides  were  able  to  specifically  block 
AR-dependent  gene  expression,  we  asked  how  various  cardiac 
glycosides  might  modulate  AR  expression.  By  RT-qPCR,  we  found 
that  these  compounds  had  little  effect  on  AR  expression,  even 
though  they  effectively  blocked  the  induction  of  the  AR-regulated 
KLK3  gene  (Fig.  5  A  and  B).  However,  by  Western  blot  analysis,  we 
found  that  the  AR  protein  was  rapidly  degraded  in  cardiac  gly¬ 
coside-treated  cells  (Fig.  5C).  Peruvoside  again  emerged  as  the 
most  potent  inducer  of  AR  degradation  with  an  estimated  IC50  of 
10-20  nM,  similar  to  its  IC50  value  based  on  cell  proliferation  (Fig. 
3 B).  AR  degradation  induced  by  cardiac  glycosides  took  place  3-6  h 
posttreatment  (Fig.  5D).  Peruvoside  was  also  most  potent  in  in¬ 
ducing  AR  degradation  in  PC3-AR  cells  (Fig.  5E). 
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Fig.  4.  Global  analysis  of  cardiac  glycosides  in  comparison 
with  AR  RNAi  on  LNCaP-abl  cells.  (A)  LNCaP-abl  cells  treated 
with  5  pM  of  individual  cardiac  glycosides  were  compared  with 
the  effect  of  AR  RNAi  (the  knockdown  efficiency  is  shown  next 
to  the  heatmap).  Significant  changes  (>twofold;  P  <  0.01)  were 
identified,  which  added  up  to  a  total  of  2,056  genes  that  were 
either  up-regulated  (red)  or  down-regulated  (blue)  on  at  least 
one  treatment  condition.  ( B )  Venn  diagrams  of  overlapped 
changes  between  AR  RNAi  and  individual  cardiac  glycosides. 
(C)  Overlapped  genes  showed  changes  largely  in  the  same 
directions. 


Cardiac  glycosides  are  well-known  inhibitors  of  the  Na+/K+ 
ATPase  via  their  direct  binding  to  the  sodium  pump  (11),  which 
was  recently  shown  to  be  sufficient  to  impair  the  induction  of  the 
IFN-p  pathway  (30).  However,  cardiac  glycoside-induced  hypoxia 
inducible  factor  (HIF)-la  degradation  in  prostate  cancer  cells 
appears  independent  of  this  mechanism  because  RNAi  knock 
down  of  the  catalytic  subunit  A  TPlal  of  the  sodium  pump  had 
little  effect  on  the  stability  of  HIF-la  (26).  We  observed  little  effect 
of  ATPlal  RNAi  in  AR  (Fig.  5F).  These  observations  suggest 
that  the  broad  anticancer  effect  of  cardiac  glycosides  may  result 
from  enhanced  degradation  of  key  cancer  gene  products.  Finally, 
we  provided  evidence  that  the  26S  proteasome  pathway  is  likely 
responsible  for  AR  degradation,  because  the  effect  could  be  fully 
suppressed  by  MG132  in  cardiac  glycoside-treated  cells  (Fig.  5G). 


Interestingly,  we  found  that  the  MG132  treatment  not  only  pre¬ 
vented  degradation  of  full-length  AR,  but  also  stabilized  a  trun¬ 
cated  AR,  which  has  been  attributed  to  alternative  splicing  of  AR 
transcripts  and  posttranslational  cleavage  of  AR  by  the  protease 
Calpaine  and  this  truncated  AR  has  been  shown  to  contribute  to 
androgen  resistance  in  prostate  cancer  cells  (31, 32).  The  induced 
proteolytical  degradation  of  AR  by  cardiac  glycosides  may  thus 
prove  to  be  an  effective  therapeutic  strategy  against  advanced 
prostate  tumors  rising  from  diverse  mechanisms. 

Discussion 

Our  present  study  elaborates  a  powerful  pathway-centric  HTS  by 
using  the  latest  deep-sequencing  technology.  This  approach  offers 
a  number  of  advantages  over  conventional  chemical  screening 
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Fig.  5.  Cardiac  glycosides  induce  rapid  AR  degrada¬ 
tion.  (A  and  B)  Effect  on  the  expression  of  KLK3  (A) 
and  AR  ( B )  in  LNCaP-abl  cells.  SDs  are  based  on  trip¬ 
licate  experiments.  (C  and  D)  Western  blot  of  AR  in 
LNCaP-abl  cells  treated  with  different  concentrations 
of  compound  for  24  h  (C)  or  with  Digoxin  for  dif¬ 
ferent  periods  of  time  (D).  (£)  Western  blot  of  AR  in 
PC3-AR  cells  treated  with  different  concentrations  of 
compound  for  24  h.  (E)  Knockdown  of  the  catalytic 
subunit  ( ATPlal )  of  the  Na+/K+  ATPase,  showing  no 
effect  on  AR  protein  in  LNCaP-abl  cells.  (G)  Pre¬ 
vention  of  induced  AR  degradation  by  the  protea¬ 
some  inhibitor  MG132.  Note  that  a  truncated  AR 
because  of  calpain-mediated  AR  cleavage  (31)  or  al¬ 
ternative  splicing  (32)  became  detectable  in  MG132- 
treated  LNCaP-abl  cells. 
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strategies.  The  approach  does  not  require  prior  identification  of 
specific  drug  targets,  thus  equally  applicable  to  both  “druggable” 
and  “nondruggable”  disease  paradigms.  This  multitarget,  path¬ 
way-centric  approach  relies  on  the  behavior  of  endogenous  genes 
(instead  of  engineered  reporter)  and  permits  identification  of  hits 
that  intervene  with  any  potential  attack  points  in  the  pathway. 
The  approach  also  overcomes  the  central  shortcoming  of  pheno¬ 
typic  screening  because  specific  gene-expression  responses  pro¬ 
vide  critical  clues  to  potential  molecular  mechanisms.  Thus,  the 
approach  described  herein  may  help  implement  a  recommended 
quick-win/fast-fail  strategy  in  early  phases  of  drug  discovery  to 
improve  the  drug  research  and  development  productivity  (4). 

The  present  screening  reidentified  a  number  of  compounds 
previously  scored  from  a  PSA  reporter  system  (16)  or  from  an  AR 
conformation  change -based  screen  (21,  22).  Interestingly,  a  group 
of  compounds  identified  from  our  screen  belongs  to  the  family  of 
cardiac  glycosides,  with  Peruvoside  showing  the  most  potent  ef¬ 
fect.  Cardiac  glycosides  have  been  previously  shown  for  their 
broad  anticancer  activities  (11).  Our  genome -wide  analysis  dem¬ 
onstrated  that  they  could  largely  mimic  AR  RNAi,  explaining 
their  antiproliferation  effects  on  both  androgen-sensitive  and  re¬ 
fractory  prostate  cancer  cells  that  still  depend  on  AR  for  growth. 

Cardiac  glycosides  have  been  best  characterized  as  inhibitors 
of  the  Na+/K+  ATPase  in  the  cell,  but  a  long  list  of  other  po¬ 
tential  mechanisms  has  also  been  documented  in  the  literature 
(11).  We  tentatively  rule  out  the  mechanism  for  induced  AR 
degradation  because  RNAi  against  the  major  subunit  of  the 
enzyme  had  little  effect  on  AR  integrity.  By  whatever  mecha¬ 
nism,  the  induction  of  AR  degradation  provides  a  plausible 
mechanism  for  the  observed  effect  of  cardiac  glycosides  in  pre¬ 
venting  prostate  cancer  among  congestive  heart  disease  patients 
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treated  with  a  widely  prescribed  cardiac  glycoside  (Dixogin) 
compared  with  untreated  groups  (12).  Specific  cardiac  glycosides 
may  therefore  be  further  developed  as  therapeutic  modalities 
against  androgen-resistant  prostate  cancer. 

Finally,  we  wish  to  emphasize  the  broad  utility  of  the  FITS2 
technology  in  both  basic  and  translational  research.  Analysis  of  a 
pathway-specific  gene  signature  coupled  with  perturbation  of  the 
pathway  by  RNAi  has  been  used  to  deduce  gene  networks  and 
crosstalk  among  Toll-like  receptors  in  response  to  diverse 
pathogens  (33).  The  FITS2  technology  would  permit  more  com¬ 
prehensive  studies  in  combination  with  genome -wide  RNAi  to 
systematically  deduce  regulatory  networks  underlying  diverse  bi¬ 
ological  pathways.  The  HTS2  technology  may  also  be  used  to  link 
SNPs  to  causal  mutations  in  human  diseases  because  the  un¬ 
derlying  RASL  assay  has  the  single  nucleotide  resolution  in 
monitoring  gene  expression  and  mRNA  isoforms  (13,  14).  Thus, 
the  HTS2  technology  offers  a  general  platform  for  large-scale 
genetics  and  chemical  genetics  studies. 

Methods 

Culture  conditions  for  various  prostate  cancer  cells,  AR  RNAi,  and  Western  blot 
analysis  were  as  described  (34,  35).  The  AR  binding  assay  was  as  previously 
described  (21).  Methods  for  chemical  screening  and  for  scoring  compound 
efficacy  and  specificity  are  detailed  in  SI  Methods. 
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